CO 

o 

CN 



CN 

< 

c3 



in 



o 



Carleman Estimate and Inverse Source Problem for Biot's Equations 
Describing Wave Propagation in Porous Media 

Mourad Bellassoued 1 2 3 and Masahiro Yamamoto 4 

1 University of Carthage, Faculty of Sciences of Bizerte, 
Department of Mathematics, 7021 Jarzouna Bizerte, Tunisia. 

2 Federation Denis-Poisson, 
FR CNRS 2964, B.P. 6759, 45067 Orleans cedex 2, France. 



<l^ ■ 3 LE STUDIUM ®, Institute for Advanced Studies, Orleans, France 



mourad.bellassoued@fsb.mu.tn 

4 Department of Mathematical Sciences, 

The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153, Japan. 

myama@ms.u-tokyo.ac.jp 

Abstract 



According to Biot's paper in 1956, by using the Lagrangian equations in classical mechanics, we 
consider a problem of the filtration of a liquid in porous elastic-deformation media whose mechanical 
behavior is described by the Lame system coupled with a hyperbolic equation. Assuming the null 
surface displacement on the whole boundary, we discuss an inverse source problem of determining a 
body force only by observation of surface traction on a suitable subdomain along a sufficiently large 
time interval. Our main result is a Holder stability estimate for the inverse source problem, which is 
proved by a new Carleman estimat for Biot's system. 



1 Introduction 



In 1956, Biot Q presented a three-dimensional theory for coupled frame-fluid wave propagation in 
fluid saturated porous media, treating the solid frame and the saturating fluid as two separate co-located 
coupled continua. Two second order coupled partial differential equations were derived from this theory. 
More precisely, let us consider an open and bounded domain S7 of R 3 with C°° boundary T = dQ, and 
let v = v(x) be the unit outward normal vector to dQ at x. Given T > 0, Biot's equation is written as: 



(1.1) 



Q U dfu s + g 12 d?uf - A^ x u s (x,t) - V (qdiviV) = F lt 

g 12 d?u s + g 2 2dfuf - V (qdiv u s ) - V (rdivu-Q = F 2 , in Q := O x (-T,T) 

with the boundary condition 

u s (x, t) = 0, u f (x, t)-u = 0, (x, t) G S := T x (-T, T) (1.2) 



and the initial condition 

(u s (x,0), <9 t u s (x,0)) = (0, 0), (V(x,0), d t uf(x,0)) =(0,0), xe!l (1.3) 



where F = (Fi, F 2 ) T is an external force with Ft = {Fj,F 2 , Ff) T , £ = 1,2, and A m a is the elliptic 
second-order linear differential operator given by 

A^ A v(a;) = / uAv(x) + (/x + A) (Vdivv(x)) 

+ (divv(x)) VA(x) + (Vv + (Vv) T ) V / u(x), xeft. (1.4) 



Throughout this paper, t and x = {x\,X2, X3) denote the time variable and the spatial variable respec- 
tively, and u s = (u{ ,u|,ti|) and u-^ = f «{,u^,w| 1 denote respectively the solid frame and fluid 
phase displacement vectors at the location x and the time t. 

Here and henceforth - T denotes the transpose of matrices under consideration. We assume that the 
Lame parameters \i and A satisfy 

H(x)>0, X(x) + fi(x) > 0, VsGO. 

The function q(x) > 0, x G O, is the dilatational coupling factor between the fluid phase and the solid 
frame. The coefficient r(x) > 0, x G $7 is the bulk modulus of the fluid phase and Q\\(x), £22 0*0 > 0, 
x G SI are the con - ected mass densities for the solid phase and the fluid phase porosity and Q12 (x) is the 
inertial coupling factor and see Horlin an Peter lTT3l . 

We assume that the sources terms are given by 

F e (x,t)=p i (x)R i (x,t), £ = 1,2, (x,t)GQ, (1.5) 

where pe G H 2 (£l) is real-valued and R? = (Rj,RJ, Rf) T satisfy 

3 

i=2 

The main subject of this paper is the inverse problem of determining p = (pi,P2) £ {H 2 {Q)) 2 uniquely 
from observed data of the displacement vector u = (u s , u* ) in a subdomain uj C fi. It is an important 
problem, for example, in mechanics to determine the source p inside a porous body from measurements 
of the slide frame and fluid phase displacements in u. 

1.1 Inverse problem 

Let u C i7 be an arbitrarily given subdomain such that dco D <9S7, i.e., w = Uflf where ^ is a 
neighborhood of F in M 3 and let R(x, t) = (R\(x, t),R2(x, t)) be appropriately given. Then we want to 
determine p(x) = (pi(x),p2(x)), x G $7, by measurements u.\ u x(-T,T)- 

Our formulation of the inverse problem requires only a finite number of observations. As for inverse 
problems for non-stationary Lame system by infinitely many boundary observations (i.e., Dirichlet-to- 
Neumann map), we refer to Rachele [36], for example. 

For the formulation with a finite number of observations, Bukhgeim and Klibanov [10] created a 
method based on a Carleman estimate and established the uniqueness for inverse problems of deter- 
mining spatially varying coefficients for scalar partial differential equations. See also Bellassoued Q], 



[2], Bellassoued and Yamamoto 0, |4), Benabdallah, Cristofol, Gaitan and Yamamoto [6], Bukhgeim 
[8 1, Bukhgeim, Cheng, Isakov and Yamamoto [9], Imanuvilov and Yamamoto lfT8l - BUI . Isakov IT241 . 
1251 . Khaidarov EH, Klibanov (291, ESI], Klibanov and Timonov t32], Klibanov and Yamamoto (331 
Yamamoto |[39l . In particular, as for inverse problems for the isotropic Lame system, we can refer to 
Ikehata, Nakamura and Yamamoto [15], Imanuvilov, Isakov and Yamamoto ED . Imanuvilov and Ya- 
mamoto [22] - [23], Isakov [24], Isakov and Kim [26]. 

A Carleman estimate is an inequality for a solution to a partial differential equation with weighted L 2 - 
norm and effectively yields the unique continuation for a partial differential equation with non-analytic 
coefficients. As a pioneering work concerning a Carleman estimate, we refer to Carleman's paper lITTI 
where what is called a Carleman estimate was proved and applied it for proving the uniqueness in the 
Cauchy problem for a two-dimensional elliptic equation. Since [11], the theory of Carleman estimates 
has been developed and we refer, for example, to Hormander lfl4l and Isakov [ 25 ] for Carleman estimates 
for functions having compact supports (that is, they and their derivatives of suitable orders vanish on the 
boundary of a domain). For Carleman estimates for functions without compact supports, we refer to Bel- 
lassoued and Yamamoto |5|, Fursikov and Imanuvilov lfT2l . lfT6ll . Lavrent'ev, Romanov and Shishat-skh 
OH, Tataru ll38l . Moreover Carleman estimates have been applied for estimating the energy and see e.g., 
Imanuvilov and Yamamoto [23], Kazemi and Klibanov 1 27], Klibanov and Malinsky [31], Klibanov and 
Timonov [32]. 

1.2 Notations and statement of main results 

In order to formulate our results, we need to introduce some notations. For xq G M 3 \r2, we define the 
following set of the scalar coefficients 

<g(m, B) = lc£C 2 (Ti), c(x) > c* > 0, x G Tl, \\c\\ c2(n) < m, Vc ' ^ ~ Xq > < 1 - # 1 , (1.7) 

where the constants m > and G (0, 1) are given. 

Assumption A. 1 

Throughout this paper, we assume that the coefficients (6ij)i <i j <2 ' A*> A, q, r G C 2 (0) satisfy the 
following conditions 

g(x) = Qu(x)q 2 2(x) - q\ 2 {x) > 0, Vie VL, 

\{x)i{x) - q 2 (x) > 0, Vx G U. (1.8) 

Let A(x) = {a.ij{x)) l<i <2 be the 2 x 2-matrix given by 




an ai2 

0-21 «22 



(1.9) 



By (1 1 -8b . we can prove that (aij(x)) 1<i <2 is a positive definite matrix on £1. 

Assumption A.2: 

Let A{x) have two distinct positive eigenvalues: ^(x), fis(x) > 0, Hi{x) / ^{x). Moreover setting, 
Ml := {q~ 1 Q22) M' we assume 

Ml,M2, V3etf(m,e). (1.10) 



Assumption A.3: 

We assume that the solution u = (u s , u/) satisfies the a priori boundedeness and regularity: 

uGtf 5 (Q), \\u\\ HHQ) < M , (1.11) 

for some positive constant Mq. 

Before stating the main result on the stability for the inverse source problem, we present Theorem 
1 . 1 on the unique existence of strong solution to (11 .It - til .21 ) with initial condition: 

u(-,0) = uo and dtu(-, 0) = ui. 

Let V(fi) = (H\n)) 3 x il(div,n), where 

#(div,fi) = {uG (L 2 (tt)f; div u G L 2 (ft)} . (1.12) 

The norm in V(£l) is chosen as follows 

ll(v\ v 2 )|| 2 m = Hv 1 !! 2 ^ + ||v 2 ||| 2(n) + ||divv 2 || 2 2(Q) , v = (v 1 , v 2 ) G v(n). 

Theorem 1.1. Let F G H l (-T,T;L 2 {ti)), (u ,Ui) G (# 2 (ft) D #o(^)) 6 x (^H^)) 6 - ^en there 
exists a unique solution u(x, t) = (u s (x, t), u^(x, £)) of ifi.iD - (fi.2D with initial data (uo, Ui) such that 

u s g c([-r,T] ; F 2 (fi) n H^(n)) n c^h^T]; f 1 ^)) n c 2 ([-r, r] ; l 2 (0)) 

u / GC 2 ([-r,T];L 2 (0)), divu / GC([-T,r];F 1 (fi))nC 1 ([-T,r];L 2 (fi)). (1.13) 
In particular there exists a constant C > such that 

\\ u \\c 2 {[-T,T];L 2 (n)) < CCII-^II^C-T.TjL 2 ^)) + II u o||h2(Q) + ||Ui ||ijl(f2))- 

Moreover, if F = 0, then the energy of the solution u = (u s , u/) given by 

E(t) = ]- I fM(a;)atu-a t u + A|divu s | 2 + 2^|e(u s )| 2 + r|divu / | 2 + 2q(divu / )(divu s ))da; 

is conserved, that is, 

E(t) = E(0), V£>0. 

Here M(x) = (^(x)/ 3 ) li < iJ < 2 and e(v) = ± (Vv + (Vv) T ). 

The proof is based on the Galerkin method and see Santos [37 ] for the case n = 2. For completeness 
we will give a proof for dimension 3 in Section 4. 

In order to formulate our stability estimates for the inverse problem we introduce some notations. 
Let d : Q — > K be the strictly convex function given by 



Set 



i?(x) = |x-x | 2 , x£(l. (1.14) 



L> 2 = maxtf(x), d 2 = mini?(x), D 2 l = D 2 -d 2 . (1.15) 



By Assumption A.2, there exist constants c* > 0, j = 1, 2, 3 such that /J,j(x) > c* > for all x G £1, 
j = 1, 2, 3. Let Cq = min{c^, C2, Cg}. We choose /3 > such that 



P + —^ry/ P < d(%, c* d z -pD z >0. (1.16) 



Here we note that since xq G" S7, such /3 > exists. 
We set 

To = * . (1.17, 

The main results of this paper can be stated as follows: 

Theorem 1.2. (Stability) Assume (A. 1), (A.2), and (A.3). Let T > T and u be the solution of dTiD-ifO) 
and (tiJD . Moreover let assume thatQj(x) := Rj(x,0) satisfy 

$j(x) ■ (x - x ) ^ forallxGTi. (1.18) 

Let M > 0. Then there exist constants C > and k G (0, 1) such that the following estimate holds: 

\\pif H ^) + \\p^ m <cunr (i.i9) 

for any p e G H 2 (Q), I = 1, 2, such that ||^||#2 < M andj^ = 0, Vy>£ = on T. Here 
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C(u) = 2^ II^ u IIh 2 (o;x(-T,T))- 

i=2 
By Theorem 1 1.21 we can readily derive the uniqueness in the inverse problem: 

Corollary 1.1. Under the assumptions in Theorem \1.2\ we have the uniqueness: 

Let u = (u s ,u-^) satisfy Biot's system rti.il )-( [01 ) such that u(x,t) = 0, (x,t) G U) x (—T,T). Then 

P\{x) = P2{x) = Ofor all x £ Q and u(x, t) = m Q. 

The remainder of the paper is organized as follows. In section 2, we give a Carleman estimate for the 
Biot's system. In section 3 we prove Theorem 1 1.21 Section 4 is devoted to the proof of Theorem ll.il 

2 Carleman estimate for Biot's system 

In this section we will prove a Carleman estimate for Biot's system, which is interesting of itself. In 
order to formulate our Carleman estimate, we introduce some notations. Let : fi — > R be the strictly 
convex function given by d 1 . 14b . where xq £ £1. 

We define two functions ip , tp : Q, x ~R — > R of class C°° by 



i{>(x,t) = \x - xof - \tf for all xG ft, -T<t<T, 



(2.1) 



ip(x,t) = e ^ ( - x ' t \ 7>0, 
where T > Tq. Therefore, by (11.17b and (11.15b . we have 

tp(x,0) > do, <p(x,±T) < do (2.2) 



with do = exp(7<i 2 ). Thus, for given r\ > 0, we can choose sufficiently small e = e{rj) such that 

<p{x, t) < d - 7] = d\ forall (x,t) G {(x,t) G Q; \t\ >T-2e}, (2.3) 
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ip(x, t) >d - -j- = d 2 for all (x, i) G {(», t) G Q; |t| < e} . 
Let (u s , u-'') satisfy Biot's system 



g u d 2 u s (x,t) + Qnd?uf - A^ x u s (x,t) - V (qdW) = Fi, 
Qi 2 d 2 u s (x,t) + g 22 d 2 u-f - V (rdiv u s ) - V (qdivu-Q = F 2 , in 



(2.4) 



The following theorem is a Carleman estimate for Biot's system (I2.4I ). 

Theorem 2.1. There exist r* > and C > such that the following estimate holds: 

f r (\V x , t u s \ 2 + |Va, t (divu s )| 2 + iV^divu')! 2 ) e^dxdt 

+ I r 3 (|u s | 2 + |divu s | 2 + (div^l 2 ) e 2rip dxdt <C I (\F\ 2 + \VF\ 2 ) e 2rv dxdt (2.5) 

for any t > t* and any solution (u s , u-f) G (H 2 (Q)) 6 to &2.4\) which is supported in a fixed compact set 
K C int(Q). 

In order to prove Theorem 12. II we use a Carleman estimate for a coupling hyperbolic system, which 
we discuss in the next subsection. 

2.1 Carleman estimate for a hyperbolic system 

First we recall the following Carleman estimate for a scalar hyperbolic equation. As for the proof, we 
refer to Bellassoued and Yamamoto [5 1, and Imanuvilov and Yamamoto GUI) for example. 

Lemma 2.1. Let c G ^(m, 6) . There exist constants C > and r* > such that the following Carleman 
estimate holds: 

C I e 2rip (t \V x , t y\ 2 + t 3 \y\ 2 )dxdt < f e 2rLp \ {d 2 t - cA)y\ 2 dxdt 

whenever y G H 2 (Q) is supported in a fixed compact set K C int(Q) and any t > t*. 
Let v = (vi,v 2 ) G (H 2 (£l)) 2 satisfy the following hyperbolic system 



d 2 v\ - 6n(x)A-ui - b 12 (x)Av 2 = g\ inQ 
d 2 v 2 - b 2 i(x)Avi - b 22 (x)Av 2 = g 2 inQ, 



(2.6) 



for g = (gi,g 2 ) G (L 2 (Q)) 2 . We assume that the matrix B(x) = {bij{x)) 1<i <2 has two distinct 
positive eigenvalues c%, c 2 G &(m, 9). Then, by Lemma 2.1, we have the following Carleman estimate. 



Lemma 2.2. There exist constants C > and r* > such that the following Carleman estimate holds 
C J e 2T v(r\V Xjt v\ 2 + T 3 \v\ 2 }dxdt< I e 2r * \g\ 2 dxdt 

for any t > t*, whenever v G H 2 (Q) is a solution of &2.6\) and supported in a fixed compact set 
K C int(Q). 

Proof. The system (12.61 ) can be written in the equivalent form 

dfv - B(x)Av = g in Q. (2.7) 

By the assumption on B{x), there exists a matrix P{x) such that 

(P- 1 BP)(x) =Diag(ci(x),c 2 (x)) = A(x), x £ tt. 

Therefore system (12.71 ) can be written in an equivalent form: 

d 2 v - A(x)Ai; = g + 3B\{x, d)v, 

where 

v(x,t) = F-^vfat), g(x,t) = P- 1 {x)g{x,t), (2.8) 

and 3B\ is a first-order differential operator. 

Since Cj G ^(m, 0) for j = 1, 2, we can apply Lemma 2.1 for the two components of v and obtain 

C I e 2rLp (t \V x , t v\ 2 + r 3 \v\ 2 )dxdt < [ e 2TLp \g\ 2 dxdt + I e 2rLp (\v\ 2 + \Vv\ 2 ) dxdt 
Jq v j Jq Jq 

and, by (I2.8I ). we easily obtain 

\v(x,t)\ < C\v(x,t)\ , \Vv(x,t)\ < C(\Vv(x,t)\ + \v(x,t)\), \g(x,t)\ < C\g(x,t)\ 
for (x, t) G Q. This completes the proof. D 

2.2 Proof of the Carleman estimate for Biot's system 

In this section, we derive a global Carleman estimate for a solutions of system (12.4b . We consider the 
6 x 6-matrix 

(Oil (x)h Qlo(x)Iz \ 
Qi2{x)h Q22{x)h ) 
Here I3 is the 3x3 identity matrix. Then by Assumption A. 1 , we have 

„,_i, v 1 / £ , 22(x)I 3 -oi 2 (x)I 3 
\ -Qi2{x)h Qu{x)h 



Let v s = divu s , v f = divu', v = (v s ,v f ) and w s = curlu s . Put G = M~ X F, F = (Fi,F 2 ) T 
and apply M _1 to system (12.4b . we obtain 

<9 t 2 u s - /xiAu s - (/ii + Ai)V(divu s )-qiV(divu^) = d + ^u s +^ divu^, inQ 

<9 t W + / u 2 Au s -r 2 V(divu-/') -q 2 V(divu s ) = G 2 + ^u s +^divu^, inQ. 

(2.10) 
Here, ^, ,£?'■, j = 0, 1 are differential operators of order j with coefficients in L°°(Q), and 

Mi = Q~ 1 H&22, Ai = ^r 1 (Ap 22 - qgi2) , <?i = 0" 1 (q£22 - r£i 2 ) 
H2 = Q" 1 VQ12, Q2 = £ -1 (q^ii - (M + A)ei 2 ), r 2 = p -1 (rgu - qgi 2 ) . (2.11) 

Henceforth &J, j = 1, ..., 4 denote some first-order operators with L°° ((^-coefficients. 
We apply div to the equations in (12- 10b . and can derive the following two equations: 

d 2 v s -anAv s -a 12 Av f = divGi + ^(V ,v s ,u s ,w s ) 

d 2 v f - a 21 Av s - a 22 Av f = divG 2 + ^ 2 {v f , v s , u s , w s ), (2.12) 

where (aij) 1<i <2 is given by (11.9I ). We apply the curl to the first equation (I2.10t to obtain 

a t 2 w s - /iiAw s = curlGi + ^ 3 (v f ,v s ,u s ,w s ) (2.13) 

and 

d t V - ftAu 8 = G X + ^ A (v f ,v s , u s , w s ). (2.14) 

Applying Lemma 12721 to system (12.121 ). we have for v = (v s ,v^) 

C f e 2rip (t \V Xtt v\ 2 + r 3 \vf\dxdt < I e 2rLp (\F\ 2 + |VF| 2 ) dxdt 

+ f e 2rip ( \u s \ 2 + |w s | 2 + |Vu s | 2 + \Vw s \ 2 ) dxdt. 



Applying Lemma |2~H to (12.131 ) and (12.141 ). we obtain 

C i e 2rv (t (V^w 8 ] 2 + r 3 |w s | 2 + r \V x , t u s \ 2 + r 3 \u s \ 2 ^dxdt 

< i e 2TV (jF| 2 + |VF| 2 ) dxdt + I e 2rip (|v| 2 + | V 
Therefore, for r sufficiently large, we obtain (I2.51 l. This completes the proof of Theorem |2~T 

3 Proof of Theorem 1.2 

In this section we prove the stability (Theorem 1.2) for the inverse source problem. 



v\ ) dxdt. 



For the proof, we apply the method in Imanuvilov and Yamamoto 11191 which modified the argument 
in iflOl and proved the stability for an inverse coefficient problem for a hyperbolic equation. For it, the 
Carleman (Theorem 2. 1) is a key. 
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3.1 Modified Carleman estimate for Biot's system 

Let cot = wx (— T, T). We modify Theorem 2. 1 for functions which vanish at ±T with first ^-derivatives. 

Lemma 3.1. There exist positive constants r#, C > and Cq > such that the following inequality 
holds: 

f r (\V X:t v s \ 2 + iV^divv 5 )! 2 + iV^divv^)! 2 ) e 2T *dxdt 

+ I r 3 (|v s | 2 + |divv s | 2 + (divv^ 2 ) e 2r<p dxdt <C f (\G\ 2 + |VG| 2 ) e 2rip dxdt 

+ Ce c ^\\A 2 H^ T) (3-D 
for any r > r* and any v = (v s , v^) G H 2 {Q) satisfying, for G = (Gi, G2) 

ft 1$ v s + ft 2 9 2 v/ - A MjA v s - V (q divv/) = Gi , 

Qi2d 2 v s + Q22d 2 vf — V (qdiv v s ) — V (rdivv^) = G2 in Q 

such that 

d{v(x, ±T) = for all x G Q, j = 0, 1. (3.3) 



Proof. Let or c u. In order to apply Carleman estimate (12.5K we introduce a cut-off function £ 
satisfying < £ < 1, £ € C°°(IR 3 ), | = 1 in fi\w° and Supp£ C fi. Let v G F 2 (Q) satisfy ((321) and 
dH. Put 

w(x,i) = £(x)v(x,t), (x,t)eQ, 

and let Qq = (fi\5J) x (— T, T). Noting that w G H 2 (Q) is compactly supported in Q and w = v in Qq 
and applying Carleman estimate (12.51 ) to w, we obtain 

/ r (iV^v*! 2 + |V a , t (divv s )| 2 + |V„ t (divv^)| 2N ) e 2r ^xcit 

+ f t 3 (|v s | 2 + |divv s | 2 + |divv^| 2 ) e 2r ^dxdt <C f (|G| 2 + |VG| 2 ) e 2rip dxdt 

+ C I \Q 2 v\ 2 e 2TV dxdt 

Jq 

for any t > t*. Here Q2 is a differential operator of order 2 whose coefficients are supported in uj. 
Therefore 

/ r (iV^Vf + |V,, t (div v s )| 2 + iV^div v^l 2 ) e 2 ^x<it 

+ / r 3 (|v s | 2 + |divv s | 2 + (divv^ 2 ) e 2r ^dxdt <C I (|G| 2 + |VG| 2 ) e 2r ^dxdt 



_i_ rv c, ° T iivii 2 

+ ^ e \\ v \\h 2 (u, t )- 



This completes the proof of the lemma. □ 



By N T;(p (v) we denote 
K, V M = f t (|V x , t v s | 2 + iV^^divv 5 )! 2 + |V^(divv^)| 2 ) e 2r ^dxdt 



I r 3 (\v s \ 2 + |divv s | 2 + |divv / | 2 ) e 2rLp dxdt (3.4) 



where v = (v s , v-^). 

Now, we recall (12.21 ) and (I2.31 l for the definition of do, rj and e and we introduce a cut-off function £ 

satisfying < C < 1, C € C°°(R) and 

C = l in(-T + 2e,T-2e), SuppC C (-T + e,T - e). (3.5) 

Finally we denote by v the function 

v(x,t) = C(t)(v s ,v^)(x,t), (x,t)eQ. (3.6) 

Lemma 3.2. There exist positive constants r*, C and Co such that the following inequality holds: 

OKA*) < / (|if + |VF| 2 ) e 2 ^dxdt + e c ^||v||^ 2(wT) + e 2d ^h\\m^ T ,T,m m 

for any r > r* and any v = (v s , v^) G (H 2 (Q)) satisfying 

Qndh s + Q 12 d 2 vf - A MjA v s - V (qdivv/) = Fi(x,t) 

Q 12 d?v s + £ 2 2^V - V (qdiv v s ) - V (rdivu') = F 2 (x,t), (x,t) G Q 

Proof. We note that v G {H 2 {Q)f and 

Pii^v* + ft2 «9 2 v/ - A M , A v s - V (qdivv/) = C(*)-Pi(ar,t) + A(v,^v), 

Pi 2 d 2 v s + ^ 22 9 t 2 v/ - V (qdiv v s ) - V (rdivv/) = ((t)F 2 (x,t) + P 2 (v, d t v), (x,t) G Q, 

where Pi and P 2 are zeroth-order operators and supported in \t\ > T — 2s. Therefore, applying Lemma 
13. ll to v and using (12.3b . we complete the proof of the lemma. □ 

3.2 Preliminary estimates 

Let ip(x, t) be the function defined by (I2.ll ). Then 

<p(x,t) = e^ x >V =: p(x)a(t), (3.7) 

where p(x) and a(t) are defined by 

p(x) = e<^ x ) > d , Vcc G and a(i) = e"^* 2 < 1, V£ G [-T,T]. (3.8) 

Next we present the following Carleman estimate of a first-order partial differential operator: 

3 

L(x,D)v = y^a,i(x)diV + ao(x)v, if!! 
i=i 
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where 

a eC(U), o = (oi,02,03)€ [C 1 ^)] 3 (3.9) 

and 

\a(x) ■ (x-x )\ > c > 0, on II (3.10) 

with a constant Co > 0. Then 

Lemma 3.3. In addition to (T3.9D and P.iOl) . we assume that ||oo|lc(m — ^ an< ^ ll^llc 1 ^) — ^, 
1 < i < 3. Then there exist constants r* > and C > such that 



r 



Kx)|V Tp(a;) dx < C / |L(x,£>MaO| 2 e 2Tp < x )ds 



for aii w e if q (O) and ah t > t*. 

The proof is direct by integration by parts and see e.g., lTT9l . 

Consider now the following system 

encfu* + Q 12 dfuf - A^ x u s (x,t) - V (qdivu*) = F x {x,t), 

Qi 2 dfvi s + Q 2 2d?uf - V (qdiv u s ) - V (rdivu-Q = F 2 (x,t), (x,t) £ Q, 



(3.11) 



with the boundary condition 

u s (x, t) = 0, u f (x, t)-v = 0, (x, i)eS (3.12) 

and the initial condition 

(u s (x,0), <9 iU s (x,0)) = (0,0), (u f (x,0), $u'M)) =(0,0), xefl, (3.13) 

where the functions Fi and F 2 are given by 

F 1 (x,t)= Pl (x)R 1 (x,t), F 2 {x,t)=p 2 {x)R 2 (x,t). (3.14) 

We introduce the following notations: 

u=(u s ,u / ), Vj(x,t) = d(u(x,t), (x,t) eQ,j = 0,1,2,3. (3.15) 

The functions Vj, j = 1, 2, 3 solve the following system 

Qndty + £ 12 d 2 vJ - A M v*(x,i) - V (qdivvj) = ^fiOM), 

2i 2 <9 2 vj + ea 2 Sf v^ - V (q div vj) - V (r div vj) = b%F 2 (x, t) , (x, t) G Q, 

with the boundary condition 

v s j(x, t) = 0, vj(x, t) ■ j/ = 0, (x, t) e S. (3.17) 

We set 
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where £(£) is given by (13.51) . We apply Lemma [3721 to obtain the following estimate: 



CMr^Vj) < / (\dlFl 2 + |Va t J F| 2 J e 2rip dxdt 

+ e ° ll v jll/P(w T ) + e \\ v j\\Hl-(-T,T;m(a.))' J=0>1; 2 >3, (3.18) 

provided that r > is large enough. 

Lemma 3.4. There exists a positive constant C > such that the following estimate 

z(x,0)\ 2 dx<C I (T\z(x,t)\ 2 + T- 1 \d t z(x,t)\ 2 )dxdt 
n Jq 

for any z £ L 2 (Q) such thatd t z G L 2 (Q). 

Proof. Let ( be the cut-off function given by (13.5I ). By direct computations, we have 

C 2 (0)\z(x,0)\ 2 dx = J° T j t ^C 2 (tMx,t)\ 2 dx^dt 

= 2 / / ( 2 (t)z(x,t)d t z(x,t)dxdt 
J-t Jn 

-2 / f ('(t)((t)\z(x,t)\ 2 dxdt. 
-t Jn 



Then we have 



z(x,0)\ 2 dx<C / (T\z(x,t)\ 2 + T' 1 \d t z(x,t)\ 2 )dxdt. 



This completes the proof of the lemma. 



□ 



Lemma 3.5. Let(f>i{x) = div (pi(x)Q e(x)) . Then there exists a constant C > such that 

2 „ 

Y, e 2 ^(\Ux)\ 2 + \VUx)\ 2 )dx 

2 „ 

< C(/V T}ip (y 2 ) + A/^(v 3 )) +Y, (\Pi\ 2 + \Vpi\ 2 ) e 2rp dx, 

provided that r is large. 

Proof. We set v^ = v| and v^ 2 ) = v^. Applying Lemma [3~4l for Zj(x,t) = e T,fi ^ x ' t ^div\'2 (x,t), 
j = 1, 2, we obtain the following inequality: 

Ct 2 [ e 2Tp V|divv^(x,0) 2 dx < r 3 [ e 2 ^ V Idiv v^(s, t) * dxdt 

Jn ~i l Jq ~i ' 

2 

+ r / e 2rv Vkdivv^^t) 2 dxcft < A/" r ,^(v 2 ). (3.19) 
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Applying Lemma l3T4l again with Wj(x, t) = e^^'^Vdiv (yV'(x, t)), we obtain 

2 2 

C [ e 2Tp y2\vdwv^(x,0) 2 dx<r I e 2rv V Ivdiv v^(x, t) 
Jn j=1 I jq j=1 I 

+ T" 1 / e 2T ^V(|Vdivd t 3 v s (x,£)| 2 + Vdivofv^(x,t 



dxdt 

J dxdi 
<AA T ^(v 2 )+AA Ti ^(v 3 ). (3.20) 



Adding (13191 ) and (15301 . we find 



i=i 



e 2Tp J] ( |divv^(x,0) + Vdivv^(x,0) ) dx < C(A/^(v 2 ) +A/^(v 3 )) . (3.21) 



Since 



\T 



M{x) (v s 2 (x,0),v f 2 (x,Q)) = (pi(x)$i(x),p 2 (x)$ 2 (x)) J , x€ft 
we have 

|^(x)| 2 + |V^(x)| 2 <C7(|v|(x,0)| 2 + |Vv|(x,0)| 2 + |v|(x,0)| 2 + |Vv|(x,0) 

+ |divv|(x,0)| 2 + |Vdivv|(a;,0)| 2 + |divv£(x,0)| 2 + |Vdivv{(x,0)| 2 

for x 6 fi. On the other hand, using (13.111 ). we obtain 



(3.22) 



vf(x,0)| 2 + |Vvf(x,0)| 2 + v£(x,0) + Vv£(x,0) < C^ (V| 2 + |V^| 2 ) , x e O. (3.23) 



Combining (13.231 ). (13.221) and (13.21b . we complete the proof of the lemma. □ 

Lemma 3.6. There exists a constant C > such that 

r f (\V Pi (x)\ 2 + \Mx)\ 2 ^ e 2 ^dx < C f (|V^(x)| 2 + \M*)\ 2 ) e 2T ^dx 

for all larger > 0, £ = 1,2. 
Proof. We have 

div((dfcp<)(x)$*(x)) = d k (t>i{x) - div (pid k $e(x)) for all fc = 1,2,3. 
Therefore 

|div ((%^)| 2 + |div ( P e$e)\ 2 ) e 2T "dx < J (|V^| 2 + |<^| 2 ) e 2r ^dx 

+ C J (\p e \ 2 + | Vp,| 2 ) e 2r ^x. (3.24) 
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Since pi = and Vpe = on the boundary F and V^ • (x — xq) / 0, we can apply Lemma 
respectively with the choice v = p£ and v = d^Pe to obtain 

r I (\d m {x)\ 2 + \p,{x)\ 2 ) e 2rp dx <cf (|div ((d kPi )$ e )\ 2 + |div (p^)| 2 ) e 2rp dx (3.25) 

for £ = 1, 2 and A; = 1,2,3. Inserting (13.241 ) into the left-hand side of (13.251 ) and choosing r > large, 
we obtain 

r J (\V Pe (x)\ 2 + \pt(x)\ 2 ) e 2rp dx <C J (|V0,(x)| 2 + \^{x)\ 2 ) e 2 ^dx. 
The proof is completed. □ 



3.3 Completion of the proof of Theorem 1.2 

By Lemmata l3.5l and l3T6l we obtain 

2 2 

T E / e2Tp(X) (l V ^)| 2 + \Plfr)f) dxKC^f (\Vpi(x)\ 2 + W{x)\ 2 ) e 2rp ^dx 

+ c(a/;^(v 2 )+a/;^(v3))- 

Therefore, choosing r > large to absorb the first term on the right-hand side into the left-hand side and 
applying (3. 19), we obtain 

T Y, I e2Tp(x) (l V ^( x )! 2 + Mx)\ 2 ) dx<C^( I (\cP t F\ 2 + |^VF| 2 ) e 2r *dxdt 

H=\ 7=2 ^ 

_L CcPvrw . II 2 I s-i 2d\T || . || 2 \ 

-r^e \\v 3 \\H*(w T )~*~ Lye HYjH/PC-t.T;// 1 ^)) J 

2 /• 

<CV/ (|V^(x)| 2 + |^(x)| 2 )e 2T ^dxdt + Ce Cor ^(u) + Ce MlT M . (3.26) 



£=1 



Then the first term of the right-hand side of (I3.261 l can be absorbed into the left-hand side if we take large 

r >0. 

Since p(x) > do, we obtain 

2 „ 

J2 / (\V P e(x)\ 2 + \pi(x)\ 2 ^j dx < Ce 2 ^- d ^ T + e c ° T ^(u) < Ce~ eT + e c ° r £ w (u). (3.27) 

At the last inequality, we used: By < d\ < do, we can choose e > such that e 2 ( d i- d o)r < e -er f or 
sufficiently large r > 0. 
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4 Well posedness of the direct problem 

This section is devoted to the study the existence, uniqueness and regularity of solutions of the following 
system: 

Ql id^M s + q 12 %vJ - A MjA u s OM) - V (qdW) = F 1 {x,t), 

£i 2 d t 2 u s + p22<9 t 2 u^-V(qdivu s )- V(rdivu') = F 2 (x,t), (x,t) € Q 

with the boundary condition 

u s (x,t) = 0, u f (x,t)-v = 0, (x, t) € S = T x (-T, T) (4.2) 

and the initial condition 

(u s (x,0), u t s (x,0)) = K, uf), (V(x,0), lifts, 0)) = (u£,u{) , xefi. (4.3) 

4.1 Function spaces 

We denote by $>{p.) the space of compactly supported, infinitely differentiable function in equipped 
with the inductive limit topology. We denote by f^'(O) the space dual to $>(£L). In general, we denote by 
X' the space dual to the function space X. We denote by (/, g) the inner product in L 2 (Q) and by (/, g) 
the value of / G X' on g € X. We use usual notations for Sobolev spaces. If X is a Banach space, then 
we denote by L p (0, T; X) the space of functions / : (0, T) — > X which are measurable, take values in 
X and satisfy: 

f Q \\f(t)\\ P X dt) = Wf\\lS{0,T;X) < 00 

for 1 < p < oo, while 

II/IIl»(o,T;X) = esssup te(0jT) ||/(t)||x < oo 

for p = oo. It is known that the space L p (0, T; X) is complete. 
We define the space 

iJ(div; n) = |u€ (L 2 (n)) 3 ; divuG L 2 (fl)| , 

equipped with the norm 

II II f\\ 11=2 , |U- ||2 V /2 

||u||H(div;n) = ^ll u llL2(n) + ll dlvu lli2 (n) J • 

Let us consider the space 

V(n) = {H 1 ^) 3 ) x #(div; ft), 

equipped with the norm 

IMIv(n) = (llu 2 !!^!^) + ||u 2 |||2 (ft) + ||divu 2 || 2 2(n) 
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4.2 Generalized solution 

We introduce the bilinear form on V(Q) by 

B(u, v) = - J ^Adiv (u s )div (v s ) + 2/j, (e(u s ) : e(v s )) + rdiv (u'jdiv (v / )) dx 

+ - q (div (u / )div (v s ) + div (v / )div (\i s )) dx (4.4) 

for any u = (u s , u-f ) G V(O), v = (v s , v-f) G V^ft). We recall that the matrix M is given by dZ9l 

Definition 4.1. We say f/iaf u = (u s ,u^) is a generalized solution of problem ( \4.1\> -< &2\) , if u G 
L 2 (0, T; V(fi)) satisfies the initial condition j[4.3\l and the following identity 

(M9 t 2 u(t),v(f)) + B(u(t),v(t)) = (F(t),v(t)) , almost all t€ (0,T) (4.5) 

for any vG £ 2 (0,T; V(J2)). 

We note that in (14.51 ) the integration is only in x. 

Lemma 4.1. Forrj > 0, we set 

B v (u,v)=B(u,v) + ?](u,v), u,v€V(n). 

Then there exists sufficiently large constant n such that the symmetric bilinear form B r] satisfies 

(i) \B v (u,v)\ < Ci||u||y (n) ||v|| v(n) , for any u,v G V(Sl), 

(ii) B v (u,u) > C 2 ||u||^^, for any u G V(Q). 

Proof. By (14.41) we obtain, for any u, v G V(fi) 

|-B(u,v)| < (||u s ||^i (n) + Hdivu^Hia^)) (||v s || H i(n) + l|divv / || L2(n) J 

<C||u||y||v||y. (4.6) 

Then for any r\, we can derive (i). 

Now, we note that for a vector u s G Hq(Q) we have the following Korn's inequality 

drt 1(n) < I e(u s ) : e(u s )dx. 
Then, for W = (div u s , div u-Q , we have 

B(u, u) > /iCi||u'||^i (n) + i /" ilW • VKdx 
where Mo is the symmetric 2 x 2-matrix given by 

Af (x) =| q | > 70I2. 

q r 
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Which implies 



B(u,u) > MCl||u s ||| 1(0) + | (||divuli 2(Q) + ||dW||f 2(n) 



1 {Q) "f" yii ulv u iiL 2 (n) ~~ yn u iiL 2 (n) 
flNl^o)- 
This completes the proof of the lemma. □ 



>C2||u||^-77||u||| 2(a) . (4.7) 



4.3 Construction of approximate solutions 

Let &/ : (L 2 (f2)) — >■ (L 2 ($7)) be the self-adjoint operator defined by 

/ A^u* + V (?dW) 
j*u =1 

\ V (gdiv u s ) + V (rdiv u / ) 

Then system (14.11 ) can be written as 

Md^u -srfu = F, (x, t) G Q (4.8) 

with initial condition 

u(x,0) = (ug(x),u£(x)), «9 t u(x,0) = K(x),u{(x)) (4.9) 

and the boundary condition 

u s (x, t)=0, u f ■ v = 0, (x, i) G E. (4.10) 

Let (wj)j>i be a sequence of solutions in (H 2 (Q) n i^Q (SI)) such that for all rn G N, wi, ...,w m are 
linearly independent and all the finite linear combinations of (wj-)j->i are dense in (i^T 2 (il)) . 
We seek approximate solutions of the problem in the form 

m 
u m{t) = ^2gjm(t)wj- (4.11) 

i=i 
The functions gj m (t) are defined by the solution of the system of ordinary differential equations 

(M3 2 u m , Wi ) + B(u m , Wi ) = {F(t), Wj) , 1 < j < m, (4.12) 

with the initial conditions 

u m (o) = u 0m ^u in (H 2 (n)nH^n))\ 

dtVLmiO) = ui m -> ui in (F 1 ^)) 6 . (4.13) 

The system (14- I2b - (I4. 13b depends on gj m (t) and therefore has a solution on some segment [0, i m ]; see 
Il35l . From a priori estimates below and the theorem on continuation of a solution we deduce that it is 
possible to take t m = T. 
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4.4 A priori estimates 

Multiplying (14.81 ) by g'j m (t) and summing over j from 1 to m, we obtain 

(Mdfu m , dtUm) +B(u m ,d t u m ) = (F(t),d t u r 
Hence 



(4.14) 



1 d 

2~dl 



r\ d 



M l l 2 d t vL m (t)\\ 2 L2m +B v {u m (t),vi m (t)) = (F(t),d t u m (t)) + ~\\n m (t)\\l 2m . (4.15) 



2dV 



»(«)• 



Let 



$ 2 (t) = ||M 1 / 2 a 4 u m (t)||2 2 +J B„(u m (t),u m (t)). 



From (I4.151 l we obtain 

\^{t)<c\\F{t)f Lm + \\d t n r . 
Integrating with respect to r from to t, we obtain 



L 2 (n) 



+ u r , 



L 2 (n) 



$ 2 (£) < C 



\ F \?m Q) + ^(0) + J (\\dtn m (r)\\ 2 L 2 {Q) + ||u m (r) 



Il 2 (Q) 



Since 



and 



®\t)>C(\\d t U m (t)\\h m + \\Urn(t"" 2 

<5> 2 (0)<C+\\u \\ 2 HHn) + ||u 1 || 2 ^ 1(n) 



V(0) 



(4.16) 

(4.17) 

(4.18) 
(4.19) 



we have from (I4.181 l 



u,, 



v(n) + ll^t u m(*)|lz,2(n) 



< R + J (\\d t u m (T)\\ 2 L2{n) + ||u^(-r)||i 2(n) ) , (4.20) 



where Rq = C + ll u o||#2(m + 1 1 u i 1 1 /mm + 11-^ Hi 2 (or ^y the Gronwall inequality, we conclude that 



u 



»»>' ittvtn) + \\dfU-m{t)\\r2(D\ < Rq 



V(0) 



L 2 (Q) 



(4.21) 



for all t e (0, T) and m > 1. 

In order to obtain the second a priori estimate, we observe that 



ll9tMo)iii a(ft) < c [i|F(o)iii 2(a) + i|u m (o)ii^ 2(ft) + i|5 tUm (o)iii 2(a) 

Indeed, multiplying (14.81) by </ m (0), summing over j and setting t = 0, we obtain 

(Md 2 u m (0),d 2 u m (0)) + J B(u m (0),d 2 u m (0)) = (F(0),d 2 u m (0)) 
Consequently, 

(Ma 2 u m (o),a 2 u m (o)) = (F(o),a 2 u m (o)) + (^u 0m ,d 2 u m (o)) 



< R\. 



(4.22) 



(4.23) 



(4.24) 
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which implies 

||3 2 u m (0)|| 2 < C (||F(0)|| 2 2(fi) + ||u 0m || 2 f2(n) ) < CR 2 . (4.25) 

Differentiating (14.141 ) with respect to t, multiplying by gj m and summing over j, we obtain the identity 

\j t [\\ Ml/2d ^m(t)\\ 2 LHn) + B v (u m (t),u m (t))] =(F(t),d?u m (t)) + ^j t \\d t u m (t)\\l 2{ny (4.26) 
Then, we conclude that 

\\d?u m (t)\\ 2 L2{Q) + \\d t u m (t)\\ 2 v <R 2 + ||d t u m (0)||| 2(n) 

+ C f (||a 2 u m (r)||2 2(n) + l^u^Cr)^) dr. (4.27) 

By (14.271) and the Gronwall inequality, we obtain 

\\d?u m (t)\\ 2 L2m + \\d t u m (t)\\ 2 v < R x . (4.28) 

Taking into consideration that u m = in X, we see 

u m e L°°(o,T-,v(n)), d t u m e L°°(o,r;y(f])), 

a 2 u m GL°°(0,r;L 2 (O)). (4.29) 

4.5 Passage to the limit 

By ( 14.291 ), we can extract a sequence from (u m ) m >o, which we denote again by (u m ) m , such that 

u m — > u in the weak-star topology in L°°(0, T; V(U)) 
dtu m — > dtu in the weak-star topology in L°°(0, T; F(O)) 
<9 2 u m -> <9 t 2 u in the weak-star topology in L°° (0, T; L 2 (0)) (4.30) 

and 

(Um, d t u m ) ->• (u, d t u) a.e., on S. 

Multiplying (14.81 ) by 6 G L 1 (0, T) and integrating, we have 

r-T f-T 



I ((Md 2 u m (t),w J )+B(u m ,w j ))6(t)dt= [ (F(t),w j )9(t)dt. (4.31) 

Jo Jo 

I B(u m ,w j )9(t)dt = - I (u m ,£/w j)9(t)dt, (4.32) 

Jo Jo 

lim / B(u m ,Wj)9(t)dt = -/ (u,j2/wj)6(t)dt = / B(u,Wj)9(t)dt. (4.33) 

l ^°°Jo Jo Jo 



'0 

On the other hand 



Thus, we obtain 



/ ((Md t 2 u(£),w,) +5(u,w j ))0(t)d*= / (F(t),Wj)6(t)dt. (4.34) 



'0 JO 
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Taking into account that Wj are dense in (i? 2 (Sl) n Hq(CI)) and therefore in V, we obtain 

(M9 f 2 u, v) + B(u(t), v(i)) = (F(t), v), i G (0, T) (4.35) 

forallveL 2 (0,T;y(^)). 

We have B(u(£), v(£)) = — (<e/u(£), v(i)) for any v G ^(^), where the application of the differential 

operator A to u is in the distributional sense in 0'(Q). Hence we obtain 

Md 2 u -£/u = F, in ^'(O), a.e. in(0, T). (4.36) 

On the other hand, <9 t 2 u, <9 t u, F G L°°(0, T; L 2 (0)). Hence (14361) holds in L°°(0, T; L 2 {Q)). 

The boundary condition (I4.21 i is satisfied by the choice of the space y(il). We prove that the initial 

conditions are satisfied. Suppose 9 G C 1 (0, T) and 6>(T) = 0. For any j we have 

J (^(u rn -u),w j \e{t)dt = -{u m {0)-u(0), Wj )6{0)- j (um^-VL^^e'Wdt. 

(4.37) 
Then, by (14.30b . we have 

lim |(u m (0)-u(0), Wi )|=0. 

m— >oo 

Since UQ m (x) = u m (0, x) and uo m — > uo, we obtain u(0) = uo, and can argue similarly for ui. 
Then, we conclude that, there exists a solution u of (14.11 ) such that 



which implies 



d t ueL°°{0,T;V{n)), and d?u G L°°(0,T; L 2 (ft)), (4.38) 



u / GC 1 (0,T;fi'(div,fi)), 9 ( V eC(0,T;L 2 ((])). (4.39) 



On the other hand 



V(qdivu s ) + Vfrdivu / J GC(0,T;L 2 (fi)), 

A„ A U S + V(qdivi/) GC(0,T;L 2 ((])). (4.40) 



Consequently, 



o 2 
A^u s GC(0,T;L 2 (fi)), A = A - — . 



Then by the elliptic regularity, u s G H^(Q.) yields 

u s eC(0,T;H 2 (n)). 
Bydivu^ G C(0,T; H 1 ^)), we see 

u s G C(0, T; # 2 (0) n Fq 1 (0)) n C^O, T; H 1 (ft)) 

i/ GC 2 (0,T;L 2 (ft)), div^ GC(0,T;^ 1 (f]))nC 1 (0,T;L 2 (O)). (4.41) 
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4.6 Uniqueness 

Let ui and 112 be two solutions to (I4.1I) - (I4.2I) with the same initial data, and set u = ui — U2. Then for 
every function v G V(Q), we have 

(Md 2 u,v) + 5(u,v) = 0, Vte(0,T). 

Since c^u G ^(fi), we may take v = dtu, and this equation can be reduced to equality 



2dt 
Then 



"M 1 / 2 ftu||| afft) + B„(u,u) 



^ " II C-AII2 

2^11^)11^(0)- 



||<9 t u(i)||| 2(n) + [|u(t)[|£ (n) < C / (||3 t u(r)|| 2 2(n) + ||u(r)||£ a(n) ) dr. 

This implies that ||u|| v"(£i) = = ||<9£u|| L 2(m and ui = U2 a.e. in Q. 
The proof of Theorem ll.ll is completed. 
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